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ABSTRACT 

The distribution of stars and gas in many galaxies is asymmetric. This so-called lopsidedness is expected to significantly affect the 
dynamics and evolution of the disc, including the star formation activity. Here, we measure the degree of lopsidedness for the gas 
distribution in a selected sample of 70 galaxies from the Westerbork H i Survey of Spiral and Irregular Galaxies. This complements 
our earlier work (Paper I) where the kinematic lopsidedness was derived for the same galaxies. The morphological lopsidedness is 
measured by performing a harmonic decomposition of the surface density maps. The amplitude of lopsidedness Ai, the fractional 
value of the first Fourier component, is typically quite high (about 0.1) within the optical disc and has a constant phase. Thus, 
lopsidedness is a common feature in galaxies and indicates a global mode. We measure Ai out to typically one to four optical radii, 
sometimes even further. This is, on average, four times larger than the distance to which lopsidedness was measured in the past using 
near-IR as a tracer for the old stellar component, and will therefore provide a new, more stringent constraint on the mechanism for 
the origin of lopsidedness. Interestingly, the value of Ai saturates beyond the optical radius. Furthermore, the plot of Ai vs. radius 
shows fluctuations which we argue are due to local spiral features. We also try to explain the physical origin of this observed disc 
lopsidedness. No clear trend is found when the degree of lopsidedness is compared to a measure of the isolation or interaction 
probability of the sample galaxies. However, this does not rule out a tidal origin if the lopsidedness is long-lived. Additionally, we 
find that the early-type galaxies tend to be more morphologically lopsided than late-type galaxies. Both results together indicate a 
tidal origin for the lopsidedness. 
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1. Introduction 

It is now known that both the stellar and the gas distribution in 
galaxies are often "lopsided". This feature was first highlighted 



Possible scenarios include tidal interactions and minor merg- 



in the classic paper by Baldwin et al. (1980i. Rix & Zaritsky 
( |1995| l analysed near-IR images of a sample of galaxies and 
characterised lopsidedness by the m = 1 mode of a Fourier anal- 
ysis. Since then, studies of larger samples of galaxies have been 
caiiied out that show that the fraction of lopsided galaxies is 
quite high and can be 50 % or more (IRichter & Sancisi||1994| 
Matthews et al. 1998 ; Bournaud et al. 2005a; Angiras et al. 2006" 
|2007j . Thus, lopsidedness is common, yet its properties are not 
as well-studied as the other common asymmetries such as spiral 
arms or bars which are characterised by the m - 2 mode. 

Lopsidedness can have a significant influence on the evolu- 
tion of the host galaxy (for a review of lopsided galaxies see e.g.. 



Jog & Combes 



2009 1. In particular, since the m = 1 mode has 



no Inner Lindblad Resonance, this mode allows the transpoit of 
matter to the central region (Blo ck et al.|r994| l. Thus, lopsided- 
ness can affect the star formation, and plays a potential role for 
fuelling the nuclear region (e.g., an AGN) ( Reichard et al. 2008 1. 

The high fraction of galaxies showing lopsidedness indi- 
cates that it is sustainable over a large period of time. However, 
its physical origin is still poorly understood. Several mecha- 
nisms might be responsible for the observed lopsided discs. 



ers ('Jog'' 1997[[Schoenmakers et al.| 1997[ [Zaritsky & Rix| 1997 



Bournaud et al.]J2005bi, asymmetric gas accretion (Bournaud 
|et al.||200 5a; Mape lli et al.||2 008 ), ram pressure from the inter- 
galactic medium (Mapelli et al.^2008j, but also an offse t of the 
stellar disc in a halo potential ( Noordermeer et al.|2001| l. 

Recent studies revealed that lopsided is coiTelated with some 
properties of the galaxies. For instance, studies of Angiras et al. 
( 2006 1 and Angiras et al. ( 2007| l indicate that the fraction of lop- 
sidedness depends on the environment. A high density group 
environment seems to favour the occurrence of lopsided discs. 
Furthermore, [Bournaud et aL] (p005a) found a coiTelation with 
galaxy type, with late-type discs being more frequently lopsided 
than early-types. In contrast, the group studies of Angiras et al.[ 
(2006 1 show that the early-type galaxies are more lopsided than 
late-types which indicates tidal encounters as a hkely mecha- 
nism for the origin of lopsidedness. 

The physical origin of lopsidedness and its life time are inter- 
related and not yet well understood. A simple kinemat ic model 
gives a lifetime of < 10^ years (Bal dwin et al.[[T980 i. A tidal 
origin which accounts for the self-gravity of the disc gives a 
somewhat longer time of 1-2 Gyr (Bournaud et al., 2005a ). A 
global mode treatment implies lopsidedness to be very long- 
lived, about the Hubble time ( |Saha et al.|2007[ ), although a non- 
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zero pattern speed would make it less long-lived Pog & Combes| 
[2009] l. 

In order to put some constraints on the physical origin of lop- 
sidedness, we here present a new study of a large sample of disc 
galaxies. Thanks to several large Hi surveys carried out in the 
last years, the number of available and useful H i maps has been 
increased enormously. We use deep medium-resolution H i maps 
of 70 galaxies from the Westerbork H i Survey of Irregular and 
Spiral Galaxies (WHISlQ. The advantage of using H i instead of 
near-IR observations is that lopsidedness can be measured out to 
much larger radii, where lopsidedness is more pronounced (Rix 
|& Z aritsky 1995 ). In fact, due to the low column density limit of 
the WHISP data, we have been able to measure the lopsidedness 
beyond the optical radius in almost all 70 galaxies. 



In a previous paper ( [van Eymeren et aL 2011 hereafter 
Paper I), the rotation curves of all sample galaxies are presented 
and the kinematic lopsidedness is analysed by comparing re- 
ceding and approaching sides. Here, we complement the results 
from this kinematic analysis by a study of morphological lopsid- 
edness at all radii. 

This paper is organised as follows: in Sect. 2, we briefly 
mention the sample selection (a more detailed description can 
be found in Paper I). In Sect. 3, the harmonic analysis is ex- 
plained. The results are presented in Sect. 4, which is followed 
by a discussion in Sect. 5. We summarise the main results and 
conclusions in Sect. 6. 



3.1. Morphological lopsidedness of the Hi gas 

The amplitude of lopsidedness in the morphology was calculated 
by decomposing the H i surface density map into Fourier com- 
ponents: 



o-(R, (f>) = aoiR) + ^ aJR)cos{m(f> - (pJR)) 



(1) 



with fl()(/?) being the mean surface density at a given radius R, 
a„, being the amplitude of the surface density harmonic coeffi- 
cient, being the azimuthal angle in the plane of the galaxy, and 
being the phase of the m''' Fourier coefficient. This harmonic 
decomposition was carried out with the GIPSY task reswri by 
using the best-fitting kinematic parameters as derived in Paper I. 
Note that all pixels along a ring are treated equally, which means 
that any small-scale structure is averaged out. The amplitudes 
were obtained up to the third Fourier component. We then cal- 
culated the normalised harmonic coefficients A,„ = a,„/a^). Ai 
represents the normalised amplitude of the m - 1 Fourier com- 
ponent denoting lopsidedness, whereas A2 represents the ampli- 
tude of an m = 2 component denoting a bar or a two-armed spiral 
feature. 

It is crucial to accurately determine the dynamic centre and 
keep it (and the systemic velocity) fixed for all rings because if 
X(), yo, and Vsys are allowed to vary, the m - 2 harmonic coef- 
ficients tend to rearrange themselves in order to minimise the 
effects of lopsidedness. In Paper I, we showed that the optical 
and dynamic centre usually agree within one beam size. 



2. Sample selection 

For the purpose of this paper we selected 70 disc galaxies from 
the WHISP survey. We applied two selection criteria: (1) the 
galaxies need to have inclinations between 20° and 75°, other- 
wise, the tilted-ring analysis will not work properly; (2) the ratio 
of the H I diameter over the beam size has to be larger than 10 in 
order to measure lopsidedness out to large radii. In order to com- 
bine reasonably high spatial resolution with sufficiently high sig- 
nal to noise, we worked on data cubes smoothed to a resolution 
of 30" X 30". As Fig. 1 in Paper I shows, the selection criteria 
do not bias our sample: the galaxies are distributed over a large 
range of B magnitudes and morphological types. They typically 
extend out to one to four optical radii, sometimes even further. 

A more detailed description of the selection process and 
some general properties of the final sample galaxies are given 
in Paper I. Details about the WHISP survey and the data reduc- 
tion process can be found in |Swaters et al.|p002[ ). 



3. Harmonic analysis 

The data analysis is based on rout ines within the Grqningen 
Image Processing System (GIPS"ypj |van der Hulst et aL]|l992| . 
A tilted-ring analysis to derive the kinematic parameters has al- 
ready been performed in Paper I. We here concentrate on de- 
scribing the harmonic decomposition, which gives us quantita- 
tive values of the morphological lopsidedness, and the calcula- 
tion of the halo perturbation. 



' http://www.astro.rug.nl/~whisp/ 

^ URL: http://www.astro.rug.nl/~gipsy/ 



3.2. Lopsidedness due to a perturbed halo 

We assume that the asymmetry arises due to the response of the 
disc to a halo distorted by a tidal encounter (Jog] 1997 1. Then, 
the perturbation potential can be obtained from the Hi surface 
density maps using the normalised amplitudes Ai, A2, and A3 
for m - 1, m - 2, and m = 3 respectively as obtained for an 
exponential disc (Jog 2000). 



In agreement with Angiras et al. ( 2006 1 we found that the 



Hi radial surface density profiles were in most cases roughly 
Gaussian in cross section. Therefore, we derived the Gaussian 
scale length for all sample galaxies: first, a radial surface density 
profile was obtained from the Hi surface density maps. We then 
fitted a Gaussian curve of exp{-{R - bf- I2R^) to the profiles 
using ?ix^ fitting technique (GIPSY task gauprof). 

The relation between ei , the perturbation parameter for the 
lopsided potential and A\ for a radial Gaussian distribution (see 
Angiras et al.|2006| l is given by 



Am 



(2) 



1 



where Ai(R) is the morphological lopsidedness and /?„ is the 
Gaussian scale length. 

The corresponding relations for m - 2 and m = 3 are given 

by 



£2 = 



A2(R) 



and 



(3) 



(4) 
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4. Results 

4.1. Morphological lopsidedness 

Figures [T] to |4] show plots of the amplitude Ai(R) and the phase 
of lopsidedness vs. R/R25, where R25 is the apparent radius 
of the optical disc. Most galaxies in our sample have H i data ex- 
tending beyond R2S (see Fig. 1, Paper I). Therefore, the Fourier 
analysis is done beyond the optical radius and in some cases 
to several times this distance. The optical disc size, R25 is typi- 
cally equal to four to five times the exponential disc scale length 
( van der Kruit & Searle|[T982) l. Thus, in the present paper, the 
Al values denoting the lopsided amplitude are measured to four 
or more disc scale lengths. This is much larger than the max- 
imum distance of 2.5 disc scale lengths to which lopsidedness 
can be studied in the near-IR ( Rix & Zaritsky|[T995 ', Bournaud' 
et al.|2005^ , where this limit is set by the high sky background. 



The previous studies showed that in the region inside of the 
optical disc, Ai increases with radius (Rix & Zaritsky 1995 



|Saha et ar]|2007[ ). This trend was also seen in the Hi studies 
of the Eridanus group galaxies that covered about four disc scale 
lengths or the optical radius ( Angiras et al. 2006). In contrast, 
in the present work we find that only about 20% of our sam- 
ple galaxies show a steady increase in the morphological lop- 
sidedness towards large radii. Good examples are UGC2455, 
UGC4173 or UGC4278. There are also examples of galaxies 
where the values of Ai(R) scatter significantly with radius (e.g., 
UGC 2953 or UGC 7989). In other galaxies, we discovered pro- 
nounced features as, e.g., in UGC 3851 or UGC 7353. And fi- 
nally, we also detected galaxies with a constant value of Ai{R) 
(e.g., UGC 3574 or UGC 7081). 

We consider /?25 to be a limiting radius between the inner 
and the outer disc. Therefore, we calculated a mean value of the 
amplitude A i (R) over all radii smaller than R25 and over all radii 
larger than R25. The results are given in Table [T] columns (4) 
and (5). The last line shows the mean of the averaged ampli- 
tudes over all sample galaxies. The mean of < Ai > for large 
radii is about 0.15, for small radii about 0.11, which means that 
the inner parts of our sample galaxies are on average slightly less 
perturbed than the outer parts. All sample galaxies showing am- 
plitudes above these values are considered as lopsided. For large 
radii, we found 30 galaxies with strong signs of lopsidedness, 
38 galaxies show values below 0.15 (for two galaxies we only 
have data within /?25)- For small radii, about 50% of all sample 
galaxies are strongly lopsided. 

Note that the typical lopsided amplitude A 1 within the optical 
disc {R < R25) of about 0.11 agrees with the value obtained by 



Bournaud et al. ( 2005a ) for the inner half of the optical disc (1.5- 



2.5 optical disc scale lengths) from near-IR studies of the OSU 
sample. Angiras et al. ( ,200 6 ) studied the Hi asymmetry in 18 
galaxies of the Eridanus group and obtained a mean larger than 
0.2 in this radial range. Therefore, our WHISP sample shows 
less lopsidedness than a group case. 

As we akeady discussed in Paper I, the existence of warps 
might affect the determination of the lopsidedness, especially 
of low to moderately inclined galaxies. However, this effect is 
not likely to be more pronounced at large radii. This is be- 
cause the amplitudes of both these features increase with radius. 
Furthermore, the amplitude of lopsidedness appears to saturate 
at radii beyond the optical disc. 

A striking result from Figs. [T] to |4] is that the phase (piiR) 
is nearly constant with radius, which indicates that m = 1 is a 
global mode. Note that 4>i(R) is given in radians. The apparently 
different values are the same angle modulo 2<p. This result will 
be discussed in detail in Sect. I 



R/R25 > 1 



R/Rjj < 1 



0.2 0.4 
<A,> 



0.2 0.4 
<A,> 



R/R25 > 1 



R/R25 < 1 



Fig. 5. A comparison of the mean values of the fractional am- 
plitudes of the first three Fourier components of the H i surface 
density distribution. Upper row: < A2 > vi. < Ai > for large and 
small radii, beyond and within the optical radius respectively. 
The < A2 > values denoting bars and two-armed spirals domi- 
nate over lopsidedness within the optical radius, for amplitudes 
> 0.2. Thus, strong bars are more likely in the inner galaxy than 
very strong lopsidedness. Lower row: < A3 > v.s. < Ai > for 
large and small radii, the < A3 > values are smaller than lop- 
sided amplitudes, especially at large radii. 



4.2. A comparison of < Ai >, < A2 >, and < A3 > 

We now want to compare the mean normalised amplitudes aver- 
aged over R/R25 > 1 and R/R25 < 1 for the first three Fourier 
components of the Hi surface density distribution. The upper 
two panels in Fig. [5] show < A2 > v^. < Ai > averaged over 
large and small radii. We see a comparable distribution for large 
radii, whereas at small radii, the m - 2 perturbations are much 
stronger. As mentioned in Sect. 



3.1 A2 denotes perturbations 



caused by a bar or spiral arms. Although spiral arms can be 
found at all radii in a disc, bars will dominate in the central parts. 
This shows that the amplitudes of the m - 2 mode are compara- 
ble with the m - \ mode at large radii, but clearly dominate at 
small radii. < A3 > is generally smaller than < Ai > independent 
of the radial range we look at (Fig. [5] lower panels). 

We calculated the mean of Ai, A2, and A3 for the two dif- 
ferent radial ranges: confirming the results from the graphs, the 
mean of A2 is lower than the mean of A\ for all radii. The mean 
of A3 is lower for all radii (see Table[T]i. 

4.3. Lopsided potential 



As shown in Sect. |3.2[ we can derive the perturbation of the 
halo up to third order from the parameters for the morphological 
lopsidedness. The values are given in the last three columns of 
Table A. 1 For completeness reasons we also list the values for 
the fractional Fourier amplitudes Ai, A2, and A3 between 1 and 
2 Gaussian scale lengths in Table |A.1| We calculated the mean 
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Fig. 1. Left panels: morphological lopsidedness Ai as a function of R/R25', right panels: phase (^i as a function of R/R25. 
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Fig. 2. Fig. [l]to be continued. 
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Fig. 3. Fig.[l]to be continued. 
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Fig. 4. Fig. [l]to be continued. 
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Table 1. Lopsidedness parameters. 
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Notes: (1) galaxy name from the UGC catalogue; (2) morphological type following the classification by |de Vaucouleurs] ( | 1 979^ ; (3) apparent 
radius from the HyperLeda database; (4) to (9) the mean values of the parameters for morphological lopsidedness as obtained in this paper; (10) 
tidal parameter as derived from Eq.[5] 



of all parameters. In agreement with the results for the morpho- 
logical lopsidedness (see Table A. 1 columns 2 to 4), the mean 
of < 62 > is higher than the mean of < e\ >, whereas the mean 
of < e3 > slightly lower than the mean of < 62 >• 

5. Discussion 

5. 1 . Variation in lopsidedness with radius 



As mentioned in Sect. |4.1[ a pronounced increase in the morpho- 
logical lopsidedness is only found in about 20% of all sample 
galaxies. This is somewhat surprising as near-IR observations at 
smaller radii show that an increase in Ai with radius seems to be 
a common feature in disc galaxies ( |Rix & Zaritsky|1995| l. Here, 
the neutral gas can be traced out to several times the apparent 
optical radius /?25- At these large radii, the increase in Ai does 
not continue, it rather appears to saturate. The origin of this dif- 
ference is not understood and presents an important clue for the 
origin of the disc lopsidedness at large radii. Unfortunately, a di- 
rect comparison of the behaviour of the stellar and gaseous com- 
ponent cannot be done here since the inner parts of the galaxies, 
on which scales Rix & Zaritsky (1995) found an increase of lop- 
sidedness, are not resolved in our H i maps. 

To get a clearer picture, we selected two galaxies and com- 
pared their global profile, the H i morphology, and the distribu- 
tions of Al and <pi with radius. UGC4173 (see Fig.|6] upper row) 
represents a type of galaxy where we found the expected be- 
haviour of increasing lopsidedness with radius. As can be seen, 
the Hi distribution is quite smooth with a clear asymmetry in 
the south-east. UGC 7256 (Fig.|6] lower row) is an example of 
a galaxy where we found strong fluctuations, so-called wiggles, 
in the distribution of Ai. The Hi intensity map shows a spiral 
galaxy with pronounced arms. The global profile has the typical 
two maxima. However, the intensities of the peaks slightly differ 
from each other, indicating lopsidedness. Furthermore, in both 
cases the phase (piiR) is nearly constant with radius indicating a 
global m - \ mode. We therefore suggest that the fluctuations in 
the distribution of A\ reflect the local Hi morphology (i.e., local 
spiral arms) superposed on top of the global lopsidedness. 

A natural way to explain the origin of these local spiral 
features would be via the swing amplification ( [Goldreich &| 
Lynden-B ell| 1 965 [ [Toomre 1 1 9 8 1 ) ) . Such fluctuations are not ap- 
parent in the near-IR studies that measure the asymmetry in the 
underlying old stellar distribution (see, e.g.. Fig. 1 in |Rix &| 
Zaritsky|[l995 ). This could be explained by the fact that a low 
dispersion component like gas has a stronger amplitude for the 
swing amplification feature than stars Pog||1992| . This is true 



even in the linear regime studied in the above swing amplifica- 
tion papers, and in a real galaxy the gas will be more likely to 
exhibit a further non-linear growth in the amplitude. 

We did a cross-check to see if this overall picture is reason- 
able. If the galaxy rotation curve shows a linear increase with 
radius, i.e., when it denotes a region of solid body rotation so 
that the angular speed is constant, the swing amplification idea 
cannot work. Hence, the net A 1 is then likely to be a pure under- 
lying global m - 1 component showing a radial increase in the 
lopsided amplitude. This is exactly confirmed for galaxies like 
UGC 4173 or UGC 4278 which show a smooth radial increase 
in Al (Fig. |2j this paper) and both have solid body rotation in 
the corresponding radial region as seen in Fig. Al, Paper I. In 
contrast, in the region of the flat rotation curve, the differential 
rotation is significant and the swing amplification is effective, re- 
sulting in the fluctuations - as confirmed for UGC 731 or UGC 
2080 which show a flat rotation (Fig. Al, Paper I) and have a 
fluctuating Ai in the same radial region (Fig.[T[ this paper). Since 
flat rotation curves are common, this explains why most of our 
sample galaxies show fluctuations in the amplitude for lopsided- 
ness Al in Hi. 

An important point to note is that while the lopsided am- 
plitude Al does show large, local fluctuations, its overall value 
increases with radius, linearly up to the optical radius and then 
somewhat slower. Thus, we can explain why lopsidedness was 
first noted in Hi in the outer parts of galaxies like M 101 
( [Baldwin et al.|1980l l. Therefore, Hi still serves as an excellent 
tracer of lopsidedness in the outer galactic discs. 

Figure |7] shows the distribution of < Ai > for small (dotted 
histogram) and large radii (solid histogram). Surprisingly, both 
distributions are comparable: most galaxies have a lopsidedness 
parameter below 0.2, only a few have higher values. Thus, while 
the lopsidedness is lower in the inner parts than beyond the op- 
tical disc (Table T}, the increase is not strong and the Ai values 
seem to saturate beyond the optical radius. This needs to be stud- 
ied by future theoretical work. 



5.2. Ttie ptiase cf>i(R) 

As already mentioned in Sect. |4.1| the phase <pi{R) is in most 
cases nearly constant with radius. This trend was already noted 
from smaller samples studied to a radial extent less than the 
optical radius ( |Rix & Zaritsky] |1995"| [Angiras et~aL] |2007| l. 
Remarkably, the phase remains nearly constant even when the 
lopsided amplitude shows a wide variation with radius as dis- 
cussed above. This indicates that m = 1 is a global mode and 



10 



J. van Eymeren et al.: Lopsidedness in WHISP galaxies: II. Morphological lopsidedness 




800 850 900 

V [km s"'] 




900 1000 1100 1200 1300 
V [km s"'] 




13"^ 30= 12" 12'^ 

R.A. (1950.0) 





5 

R/R,, 



R/R^, 



Fig. 6. Global profile (left panels). Hi intensity distribution and the distribution of Ai (middle panels), and the distribution of the 
phase (pi (right panels) of two example galaxies, UGC4173 (upper row) and UGC 7256 (lower row). 
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Fig. 7. The distribution of < A] > for small (dotted histogram) 
and large radii (solid histogram). 



denotes a perturbation in the underlying mass distribution. Since 
the stars typically contribute less than the gas to the surface den- 
sity beyond the optical radius or about four to five disc scale 
lengths as in our Galaxy (Narayan & Jog 2002), the constant 
phase most likely denotes the origin of the lopsided potential to 
be the distorted dark matter halo, to which the H i in the outer 
disc responds. This confirms our assumption made in deriving 
the perturbation potential (Sect. 3.2 1. It also makes it harder to 
support gas accretion as a mechanism for lopsidedness which 
was proposed by |Bournaud et al.] ( |2005a| l. 



5.3. Morphological vs. kinematic lopsidedness 

It has been argued that kinematic and morphological asymme- 
tries are causally connected (Jog 2002 1. However, [Komreich 
|et al.| ( |2002[ ) showed that both asymmetries are not always cor- 
related. We therefore want to find out how kinematic and mor- 
phological lopsidedness are correlated in the WHISP sample. 

A comparison of e^n and < ei >, the perturbation param- 
eter in the lopsided potential obtained by analysing the kine- 
matic and morphological data respectively, is presented in Fig. [8] 
Following the classification in Paper I, we distinguish between 
five kinematic classes: 

- Type 1: receding and approaching sides agree on all scales 
(solid triangles), 

- Type 2: constant offset of receding and approaching sides 
(solid squares), 

- Type 3: differences between receding and approaching sides 
only at large radii (open triangles), 

- Type 4: differences between receding and approaching sides 
only at small radii (open squares), 

- Type 5: receding and approaching sides change sides (stars). 

In general, it can be seen that the perturbation parameters as 
calculated from the morphological and kinematic data are com- 
parable. There are a few discrepancies where galaxies show a 
much higher value for eyn than for ej (UGC 1249, UGC 2034, 
UGC 2080, UGC 4458, UGC 4543). These galaxies are all clas- 
sified as Type 2. In case of UGC 2034, the rotation velocity could 
not be well defined as the rotation curve shows a second rise at 
higher radii. We took an average value of the first plateau as 
a lower limit for Vrot, which results in an upper limit for eidi,. 
In case of UGC 2080, we did not find any local asymmetries 
in the velocity field, but a global offset of the iso-velocity con- 
tours (see Paper I, Sect. 4.1). This results in a high e^,,,, while 
the galaxy is only moderately lopsided morphologically. Thus, 
the way the lopsidedness is measured may produce a spurious 
discrepancy between kinematic and morphological lopsidedness 
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Fig. 8. The perturbation parameter in the lopsided potential ob- 
tained by the analysis of the kinematic data < > vs. the 
one obtained from the morphological data < ei >. The sym- 
bols represent the five kinematic classes as defined in Paper I. 
Type 1: receding/approaching sides agree on all scales (solid tri- 
angles); Type 2: constant offset of receding/approaching sides 
(solid squares); Type 3: differences only at large radii (open tri- 
angles); Type 4: differences only at small radii (open squares); 
Type 5: receding and approaching sides change sides (stars). 



in some galaxies. The same seems to be true for the other three 
galaxies. 

There are also a few cases where the the perturbation param- 
eter obtained from the morphological data is higher than the one 
obtained from the kinematic data. A good example is UGC 8863 
(Paper I, Sect. 4.1) where the perturbation parameter obtained 
from the kinematic data is feeble as the discrepancies of reced- 
ing and approaching sides are small in the outer parts where the 
maximum rotation velocities have been measured. The perturba- 
tion parameter obtained from the morphological data, however, 
is much higher. 

We can conclude that morphological and kinematic lopsid- 
edness are generally comparable, which is in good agreement 
with Jog (2002) . 

5.4. Tidal origin of lopsidedness 

As mentioned in the introduction, the physical origin of lop- 
sidedness and its life time are closely connected, i.e., different 
mechanisms lead to different life times of lopsidedness. We here 
study the correlations of Ai with galaxy type and the strength 
of tidal interaction. As A; has been measured to very large dis- 
tances, we will be able to better constrain the mechanism of the 
origin of lopsidedness in comparison to earlier studies. 

In a first step, we compare the morphological lopsidedness 
< Al > averaged over large and small radii (Fig.|9] left and right 
panels respectively) to the morphological type. The early-type 
spiral galaxies seem to be more strongly lopsided, especially in 
the outer radial regions. Since a tidal interaction drives a galaxy 
towards early-type, this result of high lopsidedness in early-type 
galaxies indicates tidal encounters as the mechanism for the ori- 
gin of lopsidedness (see also Angiras et al. 2006). 

As a next step, we look at the environment of our sample 
galaxies. The degree of lopsidedness depends on the distance 
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Fig. 9. Morphological lopsidedness vs. morphological type for 
large (left panel) and small (right panel) radii. 



between the main galaxy and its companions, as well as on the 
mass ratio of the main galaxy and its companions. We therefore 
use the so-called tidal parameter to quantify the effect of tidal 
forces ( Boumaud et al.|2005a I, which is given by 



log 



Mo 



3\ 



(5) 



where M,- is the mass of each companion, Mq is the mass of the 
main galaxy, TJq is the scale length and D, is the distance of each 
companion from the main galaxy. The ratio of the masses is esti- 
mated from the ratio of absolute blue magnitudes. We take R25 /4 
as an estimate for the exponential disc scale length (van der Kruit 
& Searle 1982). The search for neighbours was done using the 



Tully Nearby Galaxies Catalogue ( |Tully & Fisher] [1988) . We 
only include companions within a radius of 2.5° and with mea- 
sured radial velocities within 500km s"^ Although almost all 
galaxies in our sample are members of groups or associations, 
many of them appear to be isolated on the basis of the above 
mentioned criteria. We find 31 galaxies to have between one and 
15 companions. The Tp values of these 31 galaxies are plotted in 
Fig. 
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vs. the parameter for morphological lopsidedness < Ai > 
averaged over large radii (large filled triangles). In a few cases, 
up to two companions do not have any data for the absolute blue 
magnitude, which gives us a lower limit for the tidal parame- 
ter (large open triangles). In two cases, only data within R25 are 
available so that we show < Ai > for small radii (small open 
triangles). 



In good agreement with |Bournaud et al. ( 2005a[ ), no corre- 
lation between < Ai > and Tp is found. Strong lopsidedness is 
as common in more isolated galaxies (smaller Tp) as in galax- 
ies with many companions. This implies that a different, internal 
mechanism for the origin of lopsidedness has to take place in 
isolated galaxies to cause them to be lopsided. Alternatively, if 
the lopsidedness is long-lived as indicated by theoretical stud- 
ies (Saha et al.|2007|l, or if it arises due to satellite accretion as 
proposed by Zaritsky & Rix| ( |1997[ | and shown by simulations 
^Boumaud et al. 2005a I, then these two factors could explain 
why isolated galaxies could be lopsided as well. If lopsidedness 
is long-lived, then there would be no correlation of Ai with a 
tidal parameter even if the origin of lopsidedness is due to a tidal 
encounter. As has been shown above, the higher lopsidedness 
seen in early-type galaxies (Fig. |9]) suggests a tidal origin. Thus, 
these two points together indicate a tidal encounter as the likely 
origin of lopsidedness in the WHISP sample studied here. 

It has to be noted that this analysis is limited by the incom- 
pleteness of the Tully Nearby Galaxies Catalogue. In some cases 
not all relevant parameters were listed in the catalogue so that 
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Fig. 10. The amplitude of morphological lopsidedness < Ai > 
averaged over large radii vs. the tidal parameter Tp (large filled 
triangles). In some cases, up to two companions could not be 
included in the calculation of Tp due to the lack of data. This 
gives a lower limit for Tp (large open triangles). In two cases, the 
H I emission could only be detected within R25 so that < A 1 > 
averaged over small radii is plotted vs. Tp (small open triangles). 
There is no clear correlation of the lopsided amplitude with the 
tidal parameter. 



these companions could not be included in the tidal parameter 
calculation. Furthermore, it is likely that not all small compan- 
ions have been detected so far. Still, this method gives a good 
indicator of the overall local environment of the galaxies, even 
though some low-mass companions might be missing. 

6. Summary 

We carried out a harmonic analysis of the Hi surface density 
distribution and measured the morphological lopsidedness in a 
sample of 70 galaxies selected from the WHISP survey. The am- 
plitude of lopsidedness, Aj, is defined to be the fractional ampli- 
tude of the first azimuthal Fourier component. The most inter- 
esting aspect of this study is that due to the high sensitivity data, 
it was possible to measure lopsidedness for the first time out to 
large radial distances. Typically, lopsidedness is measured out to 
one to four optical radii and in a very few cases even up to ten 
optical radii (Fig. 1, Paper I). 

The main results from our study can be summarised as fol- 
lows: 

- The mean amplitude of lopsidedness, < Al > within the opti- 
cal disc is high, about 0.1. The values of Ai increase radially 
within this radial range. Beyond this, the values seem to sat- 
urate. 

- The phase is remarkably constant with radius indicating that 
the lopsidedness is a global mode, as was argued earlier 
( |Sahaetal.|2007| l. 

- In a few cases, the values of A; increase with radius, but in 
most cases they show small-scale fluctuations. We argue that 
the latter are due to local spiral features arising via swing 
amplification. 

- As shown in Paper I, the galaxies can be divided into five 
groups based on their rotation velocity patterns. The most 
common types are galaxies whose velocity fields are glob- 
ally distorted (types 2 and 5). In general, kinematic and mor- 
phological lopsidedness agree quite well. 



- Early-type galaxies seem to be more lopsided than late-type 
galaxies. We did not find a correlation between lopsidedness 
and the tidal parameter indicating the strength of interactions 
for the sample galaxies. 

The results from both the morphological and kinematic studies 
indicate that tidal encounters are the main generating mechanism 
for the halo lopsidedness, to which the disc responds. 
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Appendix A: Lopsided potential 

We here present the perturbation parameters in the lopsided po- 
tential as obtained from the morphological data up to the third 
Fourier component. For a better understanding we also add the 
Gaussian scale lengths, the parameters for morphological lop- 
sidedness up to third order and the perturbation parameter in the 
lopsided potential obtained from the kinematic data in Paper I. 
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Table A.l. Lopsidedness parameters. 
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UGC /f„ < Al > < A2 > < A3 > < ekin > < ei > <€2> < 63 > 

[kpc] (l<i?„<2) (l</?„<2) (l<i?„<2) (l</?„<2) (!</?,,< 2) (1<R„<2) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


625 


16.35 


0.106 


0.157 


0.107 


0.065 


0.029 


0.052 


0.058 


731 


3.38 


0.051 


0.060 


0.073 


0.001 


0.029 


0.020 


0.043 


1249 


3.77 


0.076 


0.151 


0.138 


0.326 


0.026 


0.054 


0.084 


1256 


4.68 


0.078 


0.084 


0.056 


0.005 


0.044 


0.030 


0.037 


1281 


2.36 


0.109 


0.588 


0.098 


0.009 


0.050 


0.192 


0.062 


1317 


20.61 


0.355 


0.139 


0.141 


0.008 


0.081 


0.047 


0.095 


1501 


2.62 


0.152 


0.291 


0.126 


0.091 


0.046 


0.100 


0.070 


1913 


6.67 


0.064 


0.238 


0.083 


0.013 


0.026 


0.070 


0.055 


2034 


4.80 


0.078 


0.214 


0.056 


0.349 


0.028 


0.064 


0.033 


2080 


14.01 


0.119 


0.134 


0.051 


0.218 


0.050 


0.042 


0.034 


2455 


2.88 


0.152 


0.260 


0.062 


0.107 


0.043 


0.069 


0.037 


2800 


9.24 


0.098 


0.106 


0.138 


0.011 


0.031 


0.026 


0.071 


2855 


6.93 


0.154 


0.123 


0.099 


0.019 


0.034 


0.047 


0.061 


2953 


17.30 


0.214 


0.309 


0.098 


0.028 


0.130 


0.138 


0.073 


3273 


5.63 


0.104 


0.149 


0.048 


0.009 


0.035 


0.047 


0.027 


3371 


5.63 


0.118 


0.064 


0.114 


0.028 


0.029 


0.017 


0.065 


3574 


7.91 


0.025 


0.061 


0.061 


0.059 


0.008 


0.021 


0.036 


3580 


9.39 


0.098 


0.087 


0.057 


0.030 


0.039 


0.022 


0.037 


3734 


6.62 


0.102 


0.118 


0.108 


0.126 


0.041 


0.034 


0.061 


3851 


2.89 


0.116 


0.301 


0.214 


0.052 


0.034 


0.081 


0.136 


4173 


8.01 


0.262 


0.432 


0.225 


0.150 


0.074 


0.128 


0.138 


4278 


3.41 


0.028 


1.218 


0.050 


0.007 


0.010 


0.391 


0.031 


4284 


7.53 


0.168 


0.134 


0.078 


0.005 


0.061 


0.041 


0.047 


4458 


21.16 


0.243 


0.303 


0.144 


0.217 


0.080 


0.079 


0.077 


4543 


14.73 


0.164 


0.307 


0.065 


0.288 


0.076 


0.103 


0.041 


4838 


15.80 


0.145 


0.100 


0.149 


0.007 


0.046 


0.030 


0.094 


5079 


10.21 


0.075 


0.162 


0.085 


0.014 


0.024 


0.062 


0.053 


5251 


9.19 


0.286 


0.612 


0.161 


0.100 


0.091 


0.201 


0.096 


5253 


15.79 


0.094 


0.029 


0.040 


0.100 


0.061 


0.013 


0.029 


5532 


21.52 


0.146 


0.056 


0.086 


0.039 


0.063 


0.018 


0.058 


5685 


10.24 


0.324 


0.119 


0.051 


0.091 


0.115 


0.040 


0.032 


5717 


16.47 


0.265 


0.155 


0.100 


0.018 


0.157 


0.047 


0.060 


5721 


4.23 


0.136 


0.089 


0.143 


0.062 


0.042 


0.033 


0.090 


5789 


7.55 


0.085 


0.067 


0.047 


0.012 


0.027 


0.025 


0.032 


5829 


3.69 


0.055 


0.058 


0.096 


0.005 


0.018 


0.021 


0.056 


5918 


3.47 


0.092 


0.309 


0.044 


0.100 


0.022 


0.096 


0.028 


5997 


10.79 


0.102 


0.092 


0.125 


0.026 


0.036 


0.025 


0.075 


6225 


8.63 


0.084 


0.100 


0.053 


0.017 


0.023 


0.031 


0.032 


6446 


6.39 


0.027 


0.067 


0.044 


0.006 


0.012 


0.026 


0.030 


6537 


8.09 


0.143 


0.077 


0.080 


0.027 


0.054 


0.027 


0.052 


6787 


12.07 


0.216 


0.126 


0.072 


0.039 


0.058 


0.042 


0.047 


6937 


13.74 


0.165 


0.080 


0.111 


0.003 


0.091 


0.032 


0.079 


7081 


6.11 


0.077 


0.042 


0.079 


0.059 


0.029 


0.015 


0.047 


7090 


5.24 


0.173 


0.426 


0.062 


0.120 


0.055 


0.146 


0.039 


7095 


8.74 


0.122 


0.148 


0.134 


0.040 


0.052 


0.054 


0.085 


7151 


1.69 


0.091 


0.113 


0.034 


0.032 


0.038 


0.041 


0.021 


7256 


10.96 


0.176 


0.288 


0.180 


0.075 


0.081 


0.093 


0.108 


7321 


3.06 


0.070 


0.6900 


0.061 


0.011 


0.027 


0.222 


0.034 


7323 


3.25 


0.131 


0.077 


0.138 


0.010 


0.036 


0.022 


0.075 


7353 


9.70 


0.225 


0.273 


0.065 


0.043 


0.055 


0.097 


0.039 


7524 


3.58 


0.102 


0.087 


0.083 


0.002 


0.033 


0.028 


0.051 


7603 


2.57 


0.196 


0.101 


0.046 


0.064 


0.078 


0.031 


0.028 


7766 


13.34 


0.141 


0.039 


0.180 


0.068 


0.061 


0.012 


0.112 


7989 


11.17 


0.151 


0.356 


0.115 


0.001 


0.043 


0.107 


0.067 


8863 


11.53 


0.107 


0.195 


0.081 


0.009 


0.042 


0.058 


0.049 


9133 


42.66 


0.138 


0.071 


0.087 


0.028 


0.045 


0.022 


0.053 


9211 


4.81 


0.087 


0.080 


0.043 


0.030 


0.031 


0.024 


0.024 


9649 


3.59 


0.151 


0.081 


0.116 


0.002 


0.057 


0.024 


0.074 


9858 


15.16 


0.380 


0.250 


0.183 


0.089 


0.132 


0.077 


0.111 
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n 1 78 


u.uj i 


u.uoy 


u.uzz 


U.U4J 


U.U J 1 


11852 


31.10 


0.128 


0.386 


0.047 


0.023 


0.041 


0.104 


0.031 


11861 


9.59 


0.034 


0.097 


0.047 


0.004 


0.011 


0.031 


0.030 


11891 


4.44 


0.203 


0.054 


0.096 


0.021 


0.072 


0.017 


0.059 


12082 


4.50 


0.079 


0.239 


0.101 


0.076 


0.033 


0.075 


0.061 


12632 


4.25 


0.068 


0.050 


0.068 


0.001 


0.024 


0.017 


0.043 


12732 


7.28 


0.106 


0.092 


0.120 


0.036 


0.043 


0.032 


0.074 


12754 


3.02 


0.106 


0.042 


0.060 


0.014 


0.046 


0.013 


0.036 


mean 




0.136 


0.191 


0.096 


0.049 


0.049 


0.062 


0.059 



Notes: (1) galaxy name from the UGC catalogue; (2) Gaussian scale length as measured in this paper; (3) to (5) the mean values of the parameters 
for morphological lopsidedness averaged between 1 and 2 Gaussian scale lengths; (6) the perturbation parameter in the lopsided potential obtained 
from the kinematic data in Paper I; (7) to (9) the mean values of the perturbation parameters obtained from the morphological data in the present 
paper. 



